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ABSTRACT 



Context. During active phases of symbiotic binaries, an optically thick medium in the form of a flared disk develops around their hot 
stars. During quiescent phases, this structure is not evident. 

Aims. We propose the formation of a flared neutral disk-like structure around the rotating white dwarf (WD) in symbiotic binaries. 
Methods. We applied the wind compression model and calculated the ionization boundaries in the compressed wind from the WD 
using the equation of photoionization equilibrium. 

Results. During active phases, the compression of the enhanced wind from the rotating WD can form a neutral disk-like zone at 
the equatorial plane, while the remainder of the sphere above/below the disk is ionized. The hydrogen column density throughout 
the neutral zone and the emission measure of the ionized fraction of the wind, calculated for the mass loss rate from the WD, 
M = 2 X 10"' M0yr"' with v„ = 2 000 km s"', are consistent with those derived from observations. During quiescent phases, the 
neutral disk-like structure cannot be created because of insufficient mass loss rate. 

Conclusions. Formation of the neutral disk-like zone at the equatorial plane is connected with the enhanced wind from the rotating 
WD, observed during active phases of symbiotic binaries. This probably represents a common origin of warm pseudophotospheres, 
indicated in the spectrum of active symbiotic binaries with a high orbital inclination. 

Key words. Stars: activity - binaries: symbiotic - stars: winds, outflows 



1. Introduction 

Symbiotic stars are interacting binary systems comprising a cool 
giant as the donor star and a compact star, mostly a white dwarf 
(WD), as the accretor This composition constrains large orbital 
periods, which are typically a few years, but can be even longer 
Principal process of the interaction is the mass transfer between 
the binary components via the stellar wind from the giant. The 
following accretion heats up the WD to ^ 10^ K and increases 
its luminosity to ~ 10^-10** Lq. Such a hot and luminous WD is 
capable of ionizing neutral wind particles of both the stars, giv- 
ing rise to strong nebular radiation. The spectrum of symbiotic 
stars thus consists of three main components of radiation, two 
stellar and one nebular. Many aspects of this general view on the 
nature of s ymbi otic stars we re ori g inally pointed by, for exam- 
ple, 'Bovarchuk ( 1967), Tutukov & Yanaelson (1976), AllenJ 
Jl9m . Seaauist e t al. (1984) . Kenvon & Webbink . (19841 and 
iNuss baumer & V ogeTT ( ll987h . 

If the processes of mass loss, accretion, and ionization are 
in mutual equilibrium, the symbiotic system releases its energy 
at approximately a constant rate and spectral energy distribu- 
tion (SED). This stage is called the quiescent phase. In most 
cases, the observed energy output of ~ lO-' - 10* Lq, is be- 
lieved to be caused by stable hydrogen nuclear burning on the 
WD surface, i.e. the accreted material burns as it is accreted 
(e.g. iPaczviiski & Rudak|[T980l: lFuiimotolll982l) . During qui- 
escent phases, the main donor of particles into the circumbi- 
nary medium is represented by the wind from the giant. As a 
result, the symbiotic nebula is usually very extended, because 
the burning WD ionizes a large fraction of the giant's wind 



dSeaquist et al. 1 1984t INussbaumer & Vogel 11198 7^. Its origin al- 
lowed the mass loss r ate from the g i ant to be estimated to a few 
xlO' ^ MfTjyr"' (e.g. iMiirset et aLl Il99ll: iMikolaiewska et al.l 
I2OOI ISkopall l2O05h . Using the He method. ISkopall (120061) 
found that the mass loss rate of the hot star during quiescent 
phases is significantly smaller, being of a few xlO"^ Mq yr"'. 

Once the equilibrium between the primary processes is dis- 
turbed, e.g. by a transient increase in the accretion rate, the sym- 
biotic system changes its radiation significantly, brightens up in 
the optical by typically 1-3 magnitudes, and shows signatures 
of enhanced mass outflow for a few months to years. This stage 
is called the active phase (e.g. iKenvon II 19861) . However, a pre- 
cise mec hanism of the outbursts i s not well understood to date. 
Recently, ISokoloski et al. I {2006) suggested that nuclear shell 
burning on the WD could be enhanced by the i nflux of mate- 
rial fro m dwarf-nova-like disk instability, while iBisikalo et al. I 
(I2OO6 ') suggested a disruption of the disk as a result of varia- 
tions in the wind velocity from the giant. The SED during ac- 
tive phases is characterized by a relatively cool, A-F stellar type 
of the continuum, a decrease (or disappearance) of highly ion- 
ized lines, and an in creas e of the nebular hydrogen emission. 
Kenvon et all (11991*) and 'Mikolaiewska & Kenvon! (IT992) in- 
terpreted these features by an extended accretion disk around 
the hot star in CI Cyg and AX Per, whose plane is inclined 
to the observer with ~ 73°. To avoid the problem of an un- 
seen boundary layer, the authors speculated that the disk ma- 
terial expands out of the midplane and thus cools the bound- 
ary layer during the visual maximum. For the 1984-86 outburst 
of Z And, iFernandez-Castro et al. 1 (ll995l) also indicated a cool 
UV continuum with a signature of the Rayleigh scattering at- 
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tenuation around Ly-a, along with a strong nebular continuum, 
whose emission measure increased by a factor of 3.4 with re- 
spect to quiescence. They interpreted these characteristics in 
terms of the ejection of an optically thick shell, which sub- 
stantially blocks the UV radiation. Finally, modelling the UV 
SED of active symbiotic binaries with a high orbital inclina- 
tion confirmed the presence of such a two-temperature spec- 
trum, where the cool component is produced by a relatively 
warm stellar source radiating at » 20 OOP K, and the ho t one 
is represented by a strong nebular radiation dSkopal 112005 1. The 
former is not capable of producing the observed nebular emis- 
sion, and thus the latter signals the presence of a hot ionizing 
source 10^ K) in the system, which is not seen directly by 
the observer. In addition, a strong Rayleigh scattering attenu- 
ation of the continuum around the Ly-a line directly indicates 
a significant number of neutral hydrogen atoms on the line of 
sight, which are not present during the quiesce nt phase. This ef - 
fect was or iginally measur ed bv Fernandez-Castro et al. I (Il995h . 
Therefore, I Skopal I (l2005h suggested that there is an edge-on 
disk around the accretor, the outer flared rim of which repre- 
sents the warm pseudophotosphere. He also suggested that the 
nebula is placed above/below the disk, since it is ionized by 
the central hot star. A further important characteristic of the ac- 
tive phases is a significantly enhanced mas s loss rate from the 
hot star (e.g. Fernandez-Castro et al. I [19951 iNussbaumer et al.1 
1995 ; Crocker et al. 2002). ModelUng the broad Ha wings, 
Skopal I (12006 ) found that the mass loss rate of the active hot star 
enhances to a few times (10^^ - 10"*)Moyr"', i.e. it exceeds 
the rate indicated during quiescent phases by a factor of ^ 10. 
The ejected material is ionized by the luminous central hot star, 
which thus enhances radiation from the symbiotic nebula (see 
Tables 3 and 4 of Skopal 2005). 

These general properties of quiescent and active phases thus 
indicate a connection between the mass loss rate from the hot 
star and the level of the activity of symbiotic b inaries. Recently, 
this connection was justified bv ISkopal et al. I (2009) . who ex- 
plained the anticorrelation between the supersoft X-ray and 
UV/optical fluxes due to the variable hot star wind during dif- 
ferent stages of activity of the symbiotic star AG Dra. 

Accordingly, in this paper we investigate the idea that the en- 
hanced wind from the hot star can be responsible for the forma- 
tion of the neutral disk-Uke zone during active phases. For this 
purpose, we applied the wind compression model, as introduced 
by Biorkman & Cassinelli ( 1993), to explore how the wind from 
the WD can be compressed towards the equatorial plane due to 
its fast rotation. In Sect. 2, we outline in principle the wind com- 
pression model we use in our application, and calculate its ion- 
ization structure. Section 3 compares the model with observa- 
tions and discusses the results. Conclusions are found in Sect. 4. 

2. Ionization structure of the compressed wind 

Rotation of the hot star with radiation-driven winds leads to 
compression of the outflowing material towards the equatorial 
regions. The trajectory of each particle in the wind is determined 
by the gravity and radiation forces, which are both central forces. 
Due to conservation of the angular momentum vector, each par- 
ticle launched at the surface of the star has to move in its own 
orbital plane, which is perpendicular to this vector. The mecha- 
nism can be described by the wind compression model, in which 
the streamlines of the gas from both hemispheres of a rotat- 
ing line-driven wind are be nt towards the equatorial plane (e.g. 
iLamers & Cassinelli 1 1 1999i and references therein). The flow 
of wind particles towards the equator enhances with rotation 



and weakens with radial distance from the star's surface. Two 
forms of equatorial density enhancements can be recognized: (i) 
If the flow from the hemispheres has a supersonic component 
of the velocity perpendicular to the equatorial plane, a shock 
zone or disk will develop. This wind-compressed d isk (W CD) 
model was elaborated by iBiorkman & Cassinelli ] (119931) . (ii) 
If the streamlines do not cross the equator (the so-called non- 
crossing trajectories) we get the wi nd-compressed zone (WCZ) 
model described by llgnace et al. I (|l996). The shocks in the 
WCD model at the upper and lower boun dary of the disk can 
produce very high tempe ratures (Biorkma n & Cassinelli I [19931 : 
iKrolik & Ravmondl ll985h . As a result, the ionization structure 
in the wind would be influenced by the presence of such a hot 
disk. In our approach, we restrict calculations for the sake of 
simplicity to models with non-crossing trajectories, i.e. we ap- 
ply the WCZ model. 

2.1. Model assumptions 

The wind in our model begins at the surface of the hot WD 
with the radius 7?wd- It rotates with a velocity vvot, with the 
rotation axis perpendicular to the orbital plane. For this paper, 
we adopted 7?wd = 0.01 /?0 and v,ot - 100 - 600km s"' (see 
Sect. 2.3.1.). The photons are emitted at a blackbody WD's 
pseudophotosphere with a temperature and a radius > 
RwD- The pseudophotosphere represents the optically thick/thin 
boundary of the wind. It simulates the hot star in the model, 
which is characterized by the fundamental parameters L*, Ri,, 
and Ti,. Thus, the wind is accelerated deep inside the photo- 
sphere. As the wind mass loss rate varies between quiescent and 
active phases, is also subject to variation. In the model, we 
assume that the WD's pseudophotosphere is spherical. 

From observations, we derive the fundamental parameters 
Lh, R^, Th of the hot component in the binary (i.e. the object 
orbiting the cool giant). These parameters can deviate signifi- 
cantly from those of the hot star (i.e. L*,7?*,r* of the WD's 
pseudophotosphere as defined above), depending on the level 
of activity and the orbital inclination. For example, the effec- 
tive radius of the hot component, R'^, is defined as the radius of 
a sphere, which has the same luminosity as the entire stellar- 
radiating disk-like object, derived by modelling the UV SED 
(Skopal 2005, Eq. fT% here). During active phases, R^ can 
be as large as a few Rq for systems with a high orbital inclina- 
tion, while during quies cent phases it is clus tered around 0. 1 Rq 
only (e.g.tSkopal 2005: ISkopal et al. 11201 1[) . When viewing the 
system more from the pole, the effective radius shrinks signif- 
icantly. An example is AG Dra, for which the mean values of 
Rf = 0.063 + 0.013 7?© and 0.033 + O.OlO Rq were derived by 
model ling the SED during active and qu iescent phases, respec- 
tively (ISkopal II2005I: ISkopal et al. ll2009l) . This suggests that we 
measure R^ » for eclipsing systems and R^ ^ for 
non-eclipsing systems. Accordingly, in the model calculations 
we adopted R-i, = 0.06 Rq for active phases. Other parameters of 
the hot star and the wind that are needed for the model calcula- 
tions are introduced in Sect. 3.1. 

Finally, we neglect the effect of the orbital motion on the 
shaping of the compressed wind formation. For a typical or- 
bital period of ~ 2 years and total mass of ~ 2 Mq of a 
symbiotic binary, the orbital velocity of its components is only 
Voib ~ 30kms"'. As Vorb is significantly smaller than the termi- 
nal velocity of the wind from the WD Voo (see Sect 2.3.2), the 
momentum flux M(voo + Vorb) ~ M(vco - Voib) for a constant M 
from the WD. Therefore, the density gradient at the front and 
back side of the WD in the direction of its orbital motion will 
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also be comparable, and the expected prolongation of the com- 
pressed wind formation will probably not be significant. As a 
result, we ignore the orbital motion of the binary and calculate 
the ionization structure in the compressed wind for a stationary 
situati on and pure hydrogen gas, similar to the ISeaquist et al. I 
d 19841) calculation of the ionization structure for the quiescent 
phase. 

2.2. Density in fhe WCZ model 

According to iBiorkman & Cassinelli 1 (119931) . in a star-centred 
spherical coordinate system (r, 6, (p), the wind compression 
model assumes an azimuthal symmetry, i.e. the density of the 
wind does not depend on the azimuthal angle <p. Therefore, the 
density is a function of the radial distance r from the star's cen- 
tre and the polar angle 6, measured from the spin axis. The wind 
density distribution follows from the mass continuity equation 
as 

M Id^ ' 



47rr2yUm'«HV'r('-) \dHQ 



(1) 



where M is the mass loss rate of the star, yu^ is the mean molec- 
ular weight, ifiu is the mass of the hydrogen atom, Vi-(r) is the 
radial compone nt of the wind velocity . We us ed the y6-law wind 
as introduced bv lLamers & Cassinelli I (Il999l) . 



Vr(r) 



bR 



WD 



(2) 



where 7?wd is the origin of the wind, /3 characterises an accel- 
eration of the wind (i.e. how steep the velocity law is), and the 
parameter b is given by 



b^l- — 



(3) 



where a is the initial velocity of the wind at its origin (se e 
Sect. 2.3.2). In comparison with lBjorkman & Cassin elli"T(ll993h . 
we assume that Va, does not depend on the polar angle 9. The ge- 
ometrical factor dfx/djuo describ es the compression of the wind 
due to rotation of the star (see Biorkma n & Cassinelli I Il993t 
llgnace et al. Ill996[). A summ arv of the m odel can also be found 
in the textbook of Lam ers & Cassinelh I (|I999|) . 

Knowing the density distribution in the compressed wind 
(Eq. ([TJ), we can determine the ionization structure using the 
equation of photoionization equilibrium. 

2.3. Ionization boundaries in the WCZ model 

In this section, we calculate the ionization boundaries in the 
wind from the hot star. The ionization boundary is defined by 
the locus of points at which ionizing photons are completely 
consumed along the path outward from the ionizing star. For the 
sake of simplicity, we restrict our calculations to a wind con- 
taini ng only hydrogen atoms and assume Case B recombination 
re.g. lOsterbrock]|l974 . 

At each point of the nebula, the photoionization equilibrium 
equation can be written as 



NHo(r,e) 



f 



4nJy 
hv 



-flvdv = Np(r, e)N,{r, e)aB(T,), (4) 



where Jy is the mean intensity of radiation, Niio(r, 0), Np{r, 6) and 
A^e('", ff) are neutral hydrogen, proton and electron densities by 



number per unit volume, respectively, Oy is the ionization cross 
section for hydrogen, h is the Planck's constant, vo is the ion- 
ization frequency for hydrogen atoms, and Q'B(7^e) is the total 
hydrogenic recombination coefficient in case B. 

The left-hand side of Eq. (|4]i gives the number of ionizations 
and the right-hand one gives the number of recombinations per 
unit volume per unit time. Their equivalence expresses the con- 
dition for the ionization equilibrium. If the absorption and ge- 
ometrical dilution are taken into account, the mean intensity of 
the stellar radiation can be written as 



47Tjy = 47TW(r)By(T^)e- 



(5) 



where nBy{Ti,) is the flux at the surface of the hot star (i.e. the 
WD's pseudophotosphere, see Sect. 2.1.), Ty is the radial optical 
depth for hydrogen, and W(r) is the geometrical dilution coeffi- 
cient, which represents the ratio of the solid angle, subtended by 
the ce ntral star at the point of observation, to 'Xn (e.g. ,Gurzadvan1 
I1997h 



W{r) = i 



1 




(6) 



where 7?* is the radius of the source of the ionizing photons. 
Since the luminosity of the hot star per unit frequency interval is 
given by 



Ly = AnRinByiT^), 



we can rewrite Eq. ^ into the form 

NH<r,0) j^—aye-^^dv = Nj,(r,6)N,{r,0)as{T, 
Jv„ hynRl 



(7) 



). (8) 



Denoting the rate of photons from the hot star capable of ioniz- 
ing hydrogen as 



j^dv^^nRl 

Jv„ hy Jv„ 



nBy(T^) 
hv 



dv 



(9) 



and integrating both sides of Eq. dHJ over the radial distance r, 
we get 

Lh r Np(r,0)N,(r,0)as(T,) 



Ln n Np(r,0)N,(r,0)ai 

2nRrk i-vrr(^ 



-dr. 



(10) 



where rg is the radius of the Hi/Hii boundary at the direction 
0. To obtain the whole ionization boundary, one needs to solve 
Eq. ([Tol l for each direction, given by the polar angle 0. Due to 
the symmetry of the wind with respect to the equatorial plane 
and considering just the non-crossing streamlines of the gas, it 
is sufficient to perform calculations from the star's pole (0 = 0) 
to its equator (0 = tt/T). As in the classical Stromgren sphere 
calculations, we assume that all atoms of hydrogen with density 
distribution A^h('", 0) are ionized at the ionization radius rg (i.e. 
Np - Ns ~ Nh) and neutral outside it (i.e. Np - ~ 0). 
That means that within the ionization boundary we can take 
A^pA^e = A^^. Nussbaumer & Vogel (1987) expressed the ion- 
ization equilibrium equation in the wind with spherical density 
distribution in the form 



4n 1/" 



(r)N,(r)as{T,ydr, 



(11) 



which can be derived from Eq. (HI if we adopt the point source 
approximation (i.e. the dilution factor W = 1/4(7?* /r)^). This 
approximation is valid for Ri,/r I. However, in our treatment. 
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the ionization boundaries can be close to the central star (espe- 
cially in the equatorial plane), which requires using Eq. ( fTOl i to 
calculate correctly the H i/H ii boundaries, rg. 

Finally, we express the radial distance r in units of /?wd, and 
define u as 

M = -^. (12) 

"WD 

Then, with the further approximation of assuming a constant 
and thus a constant throughout the ionized nebula, Eq. ( fTOl i 
with the aid of Eq. ([TJ can be expressed as 



X 



du 



where ug - rgjR^Y) and the parameter X is given by 



X 



The radius of the hot star is determined in terms of its lumi- 
nosity and temperature as 



where cr is the Stefan-Boltzmann constant. Solutions of Eq. (fTST i 
for Ug at directions 9 define the H i/H ii boundary. First, we intro- 
duce some critical parameters of the rotating WD and its wind. 

2.3.1 . Rotational velocity of the WD 

Rotational velocities of the accretors in symbiotic binaries are 
not commonly known. It was observationall y found that isolated 
WDs rotate with Vro t sin / < 40kms ' (e.g. iHeber etal.lll997t 
iKoester et al. II1998I) . Therefore, it is reasonable to assume that 
accreting WDs will rotate faster than isolated WDs because of 
the angular momentum transfer from the accreting material. 

Direct evidence for a rapid rotation of accreting WDs is com- 
plicated due to a high-density circumstellar medium that shades 
photospheric absorption lines, which otherwise might be used to 
derive a possible rotation speed. An exceptio n is the case of the 
WD in the dwarf nova VW Hyi, for which ISion et a l. I d 19951) 
derived a rotational velocity of Vrotsin; ^ 600 km s"' by fitting 
the broad, shallow Si iv absorption doublet in its high-resolution 
Hubble Space Telescope spectrum. In cases where the energy of 
the accretor comes solely from the accretion process, the frac- 
tion produced by the boundary layer depends on the angular 
velocity of the accretor It was found that the observed emis- 
sion from a boundary layer in cataclysmic variables is signifi- 
cantly lower (around one-quarter) than the flux emitted by the 
accretion disk, whic h can be explained by a rapid rotat ion of 
flie WD accretor (e g. lHoare & Drew II 199 it iBelloni et all 19911: 



IVrtilek et al.lll994. Cons i dering the spinning up the accreting 



star, IPopham & Naravan I (1 19951) derived a theoretical expres- 
sion for the boundary layer luminosity. Comparing the boundary 
layer luminosity, inferred from observations, with the theoretical 
one, ISkopal et all (l2005l) found that the WD in the symbiotic 
star EG And rotates at 0.3-0.5 of its critical velocity. Further 
support for the rapid rotation of accreting WDs is provided by 
coherent oscillations in the [/-light curve of the symbiotic star 
CH Cyg with the period of 330 - 500 s, obs erved independently 
during three nights of its quiescent phase (iMikolaiewski et alTI 



119901) . The authors interpreted this periodic variability as due to 
rotation of a magnetic WD in the system. 

Accordingly, we first investigated cases with very fast ro- 
tational velocities of a WD. However, we found that v,ot ^ 
600 km s ' leads to a high compression of the wind, whose par- 
ticles can cross the equatorial plane for appropriate parame- 
ters of the WD wind. This, however, does not satisfy assump- 
tions of the WCZ model. Therefore, to calculate the ioniza- 
tion boundaries within the WCZ model, we examined cases for 
Vy„, < 600kms-i. 



(13) 2.3.2. Parameters of the hot star wind 

Modelling the broad Ha wings observed in symbiotic binaries 
by an optically thin bipolar wind from their hot components, 
ISkopal (2006) derived mass loss rates of a few xlO"*^ Mq yr"' 
and of a few x(10"^ - 10"^)Moyr"' during quiescent and 
active phases, respectively. The acceleration parameter = 

(14) 1.7-1.8 was in effect for both stages, while the wind termi- 
nal velocity during quiescence, Voo ~ 1 000 km s"' increased to 
~ 2 000km s"', during activity. 

For the initial velocity a of the wind in the model (Eq. (3)), 
we adopted the speed of sound at the interior parts of the ion- 
ized wind. According to modelling the extensive wings of emis- 

(15) sion lines from highly ionized elements by electron scattering, 
the electron temperature at the vicinity of the hot star, Te = 
30 000 - 40 000 K (see Fig. 1 of lSekeras & Skot^al 11201 2al) . cor- 
responds to the speed of sound a ~ 22kms"' (monoatomic gas 
from H and the adiabatic index of 5/3 were assumed). This value 
probably represents a lower limit, because within the acceler- 
ation zone, i.e. in the vicinity of the burning WD during outbursts 
{Ti, ^10^ K), can be considerably higher. Furthermore, because 
of the increase of the hot star luminosity, the terminal velocity 
of the wind moves from Voo ~ 1 000 to ~ 2 000 km s"' during ac- 
tive phases. Thus one can also expect an equivalent increase in 
the initial velocity of the radiatively driven wind. Accordingly, 
we adopted a - 20, 50 and 100 km s ' to calculate the Hi/Hii 
boundaries in the wind within the WCZ model. 



2.4. Results 

We calculated the H i/H ii boundaries, ug, according to Eq. ( fT3] ) 
for R^/RwD = 6 (Sect. 2.1.), v,-ot ^ 600 kms"' (Sect. 2.3.1.), 
Vco = 2 000kms-i, /? = 1.7, and a = 20,50, 100km s"' 
(Sect. 2.3.2.) with the aim of finding a range of the parameter 
X, for which a neutral disk-like zone can be created. The results 
are shown in Fig.[T| which demonstrates that for 



0.002 ^X^ 0.05 



(16) 



a neutral disk-like zone can be created at the equator of the hot 
star, flaring away from it. For X g 0.002, the ionization boundary 
tends to enclose the star, while X ^ 0.002 moves the boundary to 
farther distances with its simultaneous narrowing, until it disap- 
pears for X ^ 0.05. For a given set of parameters, different com- 
pression of the outflowing material is caused by a different WD 
rotation. Higher rotational velocities lead to a higher compres- 
sion of the wind towards the equatorial plane and to wider range 
of the parameter X, resulting in creation of the neutral disk-like 
zone (see Fig.[TJ. 

Shaping of the ionization boundaries is also sensitive to the 
wind parameters a, Vm, and /3. For a higher initial velocity a or 
terminal velocity Voo, the corresponding ionization boundary lies 
farther from the hot star. This is because higher velocities lead to 
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distance [R^d] distance [R-wp] distance [R^d] 



Fig. 1. Hi/Hii boundaries in the wind compression model (Sect. 2) at the vicinity of the hot star (soUd thick circle). They were 
calculated according to Eq. ( fT3] ) for a - 20, 50, and 100 km s ' , displayed at the top, middle, and bottom row of panels, respectively, 
with fixed v^o = 2 000km s"' and jS = 1.7. The rotational velocities, Vrot, are in kms"^ The boundaries are labelled by the corre- 
sponding parameter X. The neutral H i zone extends from the boundary in directions outside the star. The remainder part of the wind 
is ionized. Boundaries are axially symmetric with respect to the spin axis of the WD (the black circle), assumed to be perpendicular 
to the orbital plane. Distances are in units of the WD radius. 



lower densities in the wind, which moves the point of the balance 
between the flux of ionizing photons and neutral atoms to farther 
distances. On the other hand, a higher value of the parameter /3 
accelerates the wind from a to Vt» more slowly, leading to higher 
densities in the wind and thus to a closer position of the neutral 
zone to the hot star, which is characterized by a larger opening 
angle. 



The range of the ionization parameter X, for which the neu- 
tral disk-like zone can be created, is important for comparison 
with observations. 

3. Comparison with observations 

Application of the wind compression model to the hot compo- 
nents of symbiotic binaries during active phases revealed the 
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-8.5 -8 -7.5 -7 -6.5 -8 -5.5 

logM[Moyr"'] 



Fig. 2. Parameter X as a function of M from the hot star, as 
given by physical parameters of the selected symbiotic stars dur- 
ing quiescent (left) and active (right, keys) phases. Their fits 
with Eq. ( [TtT i are denoted by the solid lines. The grey belt cor- 
responds lo X - 0.002 - 0.05, which leads to creation of the 
neutral zone around the hot star in the form of a flared disk (see 
Fig.[I] Sects. 2.3. and 2.4.). 



possibility of forming a neutral disk-like zone around the hot 
star. This confirmed theoretically the previous suggestion based 
on modelling the UV/near-IR SED that there is a neutral flared 
disk surround ing the active hot star at the orbital plane (see 
Fig. 27 of Skopal 2005). The model is also consistent with the 
ionization structure derived from the photometric and spectro- 
scopic analysis of the eclipsing symbiotic binary AX Pe r during 
its recent active phase (see Fig. 9 of iSkopal et al. II20TTI) . 

In the following sections, we verify the model by compar- 
ing its basic properties with those obtained independently from 
observations. In Sect. 3.1., we investigate the dependence of the 
wind compression, given by the parameter X, on the mass loss 
rate from the hot star. The aim is to determine if the presence 
or absence of the neutral disk around the hot stars during active 
or quiescent phases is consistent with the model predictions. In 
Sect. 3.2., we compare the model column density of the neutral 
hydrogen, calculated throughout the neutral zone, with that de- 
rived from observations. Finally, in Sect. 3.3., we compare the 
observed emission measure to that given by the model. 

3.1. Wind compression as a function of M 

According to Eq. (fT3T l. the ionization structure in WCZ model 
is given by the parameter X. It reflects a degree of compression 
of the wind (Fig. [T]i- The changes of the X parameter are dom- 
inated by significant changes in the mass loss rate from the hot 
star during different phases of the star's activity (Sect. 2.3.2.). 
Therefore, to probe if the presence or absence of the neutral 
disk in a symbiotic system is consistent with the WCZ model, 
we investigate a dependence of X on M for physical parameters 
of selected symbiotic stars during quiescent and active phases. 
According to Eq. ( fT4l i. the value of the parameter X is given by 
the hot star parameters Lh, R* and those of its wind, M, Voo. 



3.1 .1 . X(M) relationship during quiescent pliases 

For a quiescent phase, we selected ten objects with available 
Lh, R'^, Th of their hot components, as derived by modelling 
the UV/near-IR SEDs (EG And, Z And, CI Cyg, LT Del, 
YY Her , RW Hya, SY Mus, AG Peg and AX Per; see Table 3 of 
'Skopal >'2b05) and X-ray/near-IR SED (AG Dra; see Table 3 of 
Skopal et al. 2009). Lh is determined by Lh an d Th. Pa r ameter s 
of the wind were taken from Table 1 of ' Skopal I ( |2006|) . 
For AG Dra and LT Del, we adopted average values Voo = 
1 000 km s- ' and M = 5 X 1 Mq yr- 1 . In the case of EG And, 
w e used Vco = 900kms ~' and M - 5 x 10"' Mq yr"', according 
to lMiirset et al. I (1 19971 ). 

To get a rough estimate of the X(M) quantities, we first used 
fundamental parameters L^, R^, as derived from observa- 
tions. The results are plotted in Fig. |2]by full grey circles. The 
'plus' sign coiTesponds to AG Dra, which is seen more from 
the pole, and thus its R^ (see Sect. 2.1.). This implies 

that the other nine systems (all with a high orbital inclination) 
have R^ > which corresponds to a lower value of X (see 
Eq. (fT4ll). To improve the results, we determined the X(M) quan- 
tities for the hot star parameters Li,,Ri,,Ti, as follows, (i) We 
put s Th = 160 000 K, as derived by modelling the X- 
ray/UV SED of AG Dra. Such a high temperature is also sup- 
ported by a direct X-ray observation of a sample of symbiotic 
stars with a pronounced supersoft X-ray (0. 1-0.4 keV) radiation 
that is si milar in th e profile with that of AG Dra (see Fig. 2 of 
iMurset et'anil997h . (ii) We determined 7?* and Lh for L* = Lh- 
Corresponding X(M) values are plotted in Fig.|2]by open circles. 

According to Eq. (fT4l l. the X(M) relationship can be written 
in the form 

log(X) = C-21og(M), (17) 

where C = log(9.6 x IO'^^/^^^/T^^L h) + 21og(vco) for 
QfB(H,20000K) = 1.43xl0-i^cm3s-' (e.g. lHummer & SeatonI 
1963) and given parameters. Fits to filled circles and to open 
circles with the 'plus' sign in Fig. |2] correspond to C = -14.81 
and -14.26, respectively (for M in MQyr"'). This shows that 
the quantities of X for the hot star parameters (L*, L*) are a 
factor of ~ 3.5 higher than those determined using the observed 
parameters (Lh, R^, Lh). Thus, regardless of the selected param- 
eters, both sets of the X(M) data are located well above the grey 
belt in the figure, i.e. without the possibility of creation of the 
flared neutral disk-like zone, as predicted by the WCZ model 
(Sect. 2.4., relation (fT6ll). In other words, the mass loss rate from 
the hot star during quiescent phases is too small to give rise to 
a neutral zone of hydrogen, because the rate of ionizing photons 
is sufficiently high to ionize even the compressed fraction of the 
wind. 

3.1.2. X{M) relationsiiip during active pliases 

Determination of the hot star parameters during active phases 
for systems with a high orbital inclination is complicated by 
the presence of the edge-on disk, which blocks the original ra- 
diation from the hot star in the direction of the observer (see 
Sect. 5.3.6. of Skopal 2005). Therefore, indirect methods based 
on the ratio of nebular line fluxes and/or the presence o f the high- 
est io nization states in the spectrum (iMiirset & NussbaumeTI 
1 19941) are employed to estimate the hot star temperature dur- 
ing active phases. For example, using the ratio of He 11 4686 A 
and H/3 line fluxes from the spectr a of the symbiotic proto type 
Z And during its 2000-03 outburst, ISokoloski et al. I (l2006h de- 
rived Lh = 92 000 + 20000 to 160000 + 35 000 K. During the 
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smaller, 1997 outburst, both methods indicated an increase in Th 
to 180000 K. 

More trustworthy estimates of Tj, can be achieved only for 
non-eclipsing systems with the aid of multiwavelength mod- 
elling of the supersoft X-ray/UV SED, which covers the total 
spectrum emitted by the burning WD. An applicable system here 
is AG Dra, the strongest supersoft source among classical sym- 
biotic stars. Although the supersoft component disappears dur- 
ing its active phases, the very high far-UV fluxes and a strong 
nebular continuum constrain Th to be > 150 - 180kK, depend - 
ing on the outburst stage (see lSkopalll2005l:ISkopal et al. Il2009h . 
As the orbital inclination of AG Dra is low (; ^ 30 - 60°, 
iMikolaiewska et al. 1995; Schmid & Schild Ill997h . we can di- 
rectly see the WD's pseudophotosphere and thus can assume that 
its T-i, X! Th and R-,, (see also Sect. 2.1). Therefore, we 

adopted the average hot star parameters, which were derived 
from the model SEDs of AG Dra, = 170000 + 7000 K, 
and R^ = 0.06 + 0.01%, (i.e. = 2 700 ± 700 L© and 
Lii = ( 1 .4 + 0.4) X 1 0"*^ s" ' ), as representative values for the active 
phases of AG Dra, Z And, AX Per, AE Ara, and AR Pav. The pa- 
rameters VjoMid M for these systems are available in Table 1 of 
ISkopall (I2006h and Table 7 of lSkopal et all (120111) . The X(M) 
quantities of active symbiotic stars lie well within the possibility 
of creation of the flared neutral disk-like zone predicted by the 
WCZ model (the grey belt in Fig.|2]l. In other words, during ac- 
tive phases the mass loss rate is so high that the rate of ionizing 
photons is not capable of ionizing the compressed wind around 
the hot star equator. 

Figure |2] suggests that the presence or absence of the neu- 
tral disk-like zone in a symbiotic system during active or quies- 
cent phases, respectively, is given by the value of the parameter 
X. During quiescent phases X x I, while during active phases 
X =s 0.01. This large difference results from the very different 
values of M in quiescent and active phases (Sect. 2.3.2.) and the 
proportionality of X to M^^. Thus, we can conclude that creation 
of the neutral disk-like structure, indicated in the spectrum of ac- 
tive symbiotic binaries, is consistent with the WCZ model, i.e. it 
can be a result of a compression of the hot star wind, which is 
significantly enhanced during active phases. 

3.2. Column density of the neutral hydrogen 

In this section, we verify the model by comparing the column 
density of the neutral hydrogen atoms, which was calculated 
throughout the neutral disk-like zone, with that derived from ob- 
servations during active phases of symbiotic binaries. 

Modelling the SED of symbiotic binaries during active 
phases showed the presence of a warm (1 - 2 x 10** K) disk-like 
pseudophotosphere (Sect. 1) and a large amount of the neutral 
hydrogen near to the orbital plane of the binary. Column densi- 
ties derived from the ultraviolet spectra are between ~ 10^' and 
afewxlO^^ cm"^ (Skopal 2005). A large amount of the neutral 
material can be identified even around the orbital phase ~ 0.5, 
when the hot component is in front of the cool giant. Figure [3] 
demonstrates this case for BF Cyg. During the active phase, the 
radiation of the warm pseudophotosphere was Rayleigh attenu- 
ated with «H = 4.8 X 10^^ cm"^ (the grey part in the figure), while 
during quiescence, the hot component radiates at Th > 10^ K 
with no signatures of the Rayleigh scattering and the iron cur- 
tain absorptions. This demonstrates that the presence/absence of 
the large amount of the neutral hydrogen on the line of sight 
is connected with intrinsic variations of the circumstellar mat- 
ter around the hot star during active/quiescent phases. A con- 
tribution from the interstellar matter (ISM) in the direction to 
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Fig. 3. Ultraviolet SED of BF Cyg during its 1990 optical max- 
imum (mid) and the following quiescent phase (bottom). Times 
of observations are denoted by arrows at the top panel with the 
light curves. In both cases, the binary was seen from the same 
position with the hot star in front (orbital phase if ~ 0.6). The 
dashed and the thin solid line represent the continuum from the 
nebula and the hot star, respectively. The thick solid line is the 
model SED. The grey surface marks the Rayleigh attenuated 
continuum with «h = 4.8 x 10^^ cm"^ , and strong blended ab- 
sorptions of the iron curtain are denoted by 'i.e.'. 



0.35 mag (iMtirset et all 



1 .7 X 10-1 cm-2 (e.g. lDiplas & Savage 119941 )? 



BF C yg can be neglected. For E^- 
n^gT), nH(ISM) 
which is difficult to measure on the low-resolution International 
Ultraviolet Explorer spectra. 

In the WCZ model, we calculated the column density of the 
neutral hydrogen atoms in the radial direction 9 throughout the 
neutral zone as a function of r. 



Ji>i„ 



(18) 



where the density distribution in the wind, Nuir, 6), is given 
by Eq. ([TJ. Figure |4] shows two riniO, r > rg) functions cor- 
responding to ionization boundaries calculated for X - 0.01, 
v'lot = 300kms ', a - 20kms"' (see the top right panel in 
Fig. [TJ, and two different orbital inclinations i = 9 = 80° and 
90°. The edge-on view (9 - 90°) gives a maximum value of 
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Table 1. Column densities «h, effective radii of the hot com- 
ponents R^, disk radii R^, and emission measures EM de- 
rived from observa tions. Data are from .Skopal . (.2005i) and 
ISkonal et al.l(l2Qllh . 



Object 


log «H 

(cm ) 


Rf 

(«o) 






EM 
cm ) 


BFCyg 


— ^ , ,-,_i_n 1 c 

22.68!«;^g 


7.1 


13.0 


24±8.3 


< 120 
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23.00!°]^ 


0.67 
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22±8.0 


5.3 


YYHer 
AR Pav 
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23 14+""* 
22.65!°]^ 
22.84!°;^^ 


0.89 
1.90 
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6.2 


I. 63 
3.48 
0.77 

II. 3 


20±8.0 
21±8.0 


7.2 
17 
3.6 
3.5 
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3.8 
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22±2.0 


8.2 






7.0 


12.8 




4.0 
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4.4 



according to Eq. il9t . ^' from the eclipse profile (see text). 



the column density, since we are looking through the densest 
parts of the compressed wind. Calculations were performed for 

= 1.7x lO^KandL* = 3 x 10^ L©, i.e. Lh = 1.6x lO'^^s-' 
(Sect. 3.1.2.), and the wind pai-ameters, M = 2 x 10"'' Mq yr"\ 
Voo = 2 000km s"', which represent mean values measured dur- 
ing active phases (see Table 1 of Skopal 2006). 

To compare the observed values of ne with the model (i.e. 
riniO, r - Ri), see Eq. (fTSTi). we need to know the radius of the 
optically thick disk, R^, whose outer rim represents the warm 
pseudophotosphere indicated during active phases of symbiotic 
binaries with a high orbital inclination (Sect. 5.3.5 and Fig. 27 in 
ISkop al 2005). Therefore, a rough estimate of R^ can be obtained 
from the effective radius, R'^ , of the hot components derived 
from the SED. According to its definition (here in Sect. 2.1.), we 
can write 



InRilHa-T^, 



(19) 



where H is the height of the flared disk (see also Fig. |6]l. For 
a can onical value of ////?d = 0.3 during active phases (Skopal 
120061) . Ri ~ 2 X Rf . Table 1 summarizes va lues oi Rf fr om 
Table 4 of lSkopall (12Q05I) and Table 8 of Skopal etaLj (1201 l b. 

The disk radius can also be estimated directly from eclipses, 
which deve lop in optical light curves duri ng act i ve ph a ses (e.g. 
Belvakina] ll99U Munari 1992; Brandi etal.l 120051: ISkopall 
2008 . Fig. 4 here). For the 1991 eclipse of BE Cyg, the radius 
of the eclipsed object in units of the binary components sepa- 
ration, A, was R^IA = 0.051 + 0.018, which yields R^{= R^) = 
24 + 8.3 Rq for the total mass of the bi nary Mt - 2.5 Mq and the 
orbital period of 757.3 days ( Skopal e t al. II 19971 and references 
therein). CI Cyg: The first two contact times of its 1975 eclipse, 
Ti^JDl 442 634.5 + 1 . 5, Tj^ JD 2 442 646.8 + 2.0 (from the 
light curve of Belvakin all 19911) an d its middle, 7Dec1.(1975.8) = 
2442691.6 ± 0.8 (iSkopallll998h . orbital period of 853.8 days 
dFekelet al.ll2000h and Mt = 2.0 Mq (Kenyon et al. 1991J give 

= 22 + 2>Rq. YY Her: The contact times of the eclips e 
that developed during the 1993- 94 active phase (ISkopal Il2005l) . 
the orbital period of 574.6 days jWiecek et al. 11201 Oh for Mt = 
2.0 Mq gives A = 361 Rq and thus R^ ^ 20 + SRq. AR Pav: 
Durin g epochs 59 - 63 ( 1987 - 1995) of its historical light 
curve. ISkopal et al. I (l2000h estimated R^/A ^ 0.05 ± 0.02 (see 
Fig. 1 there), which for A = 419 (ISchild et all 1200 ll) im- 



24.0 
23.5 
c." 23.0 

's 

o 

^ 22.5 
fi 

60 

.2 22.0 
21.5 
21.0 



model: X = 0.01, v^ot = 300 km s"\ a = 20 km s"^ 




YY 

CI Cyg O 

AX Per □ 

BF Cyg ■ 

AR Pav A 

10 20 30 

radial distance [Rq] 



9 = 80° 



40 



50 



Fig. 4. Comparison of the model hydrogen column densities 
nH(90°,r > rg) and «H(80°,r > rg) (solid lines), calculated ac- 
cording to Eq. ( ITSl l. with those derived from observations of se- 
lected eclipsing objects, plotted at the distance r - R^ from the 
WD centre (Table 1). 



plies «d = 21 + 8%. AX Per: Recentlv. ISkopal et al. I (1201 ll) 
derived R^ - 21 + 2 Rq from the 2009 eclipse, which developed 
during the 2007-10 active phase. The eclipse profile was asym- 
metric, with the first two contact times corresponding to R^ - 
22.3+2 Rq, in particular agreement with the disk radius derived 
from the model SED {R^^ - IIRq and Eq. ([T9]l). It is important 
to note that the radius of the warm disk-like shell, derived from 
the eclipse profile, represents an upper limit because of a mea- 
surable contribution of the ionized part of the hot star wind at 
radial distances of ~ 20% (EM{r > 20 Rq) ~ 10" cm"^^ in the 
model). This also causes large uncertainties of /?d from eclipses. 

Figure |4| compares quantities of «h derived from observa- 
tions of symbiotic stars with a high orbital inclination, BF Cyg, 
CI Cyg, YY Her, AR Pav, and AX Per. We used column densities 
derived from the Rayleigh scattering of the far-UV continuum 
photons around the Ly-ff line on the neutral atoms of hydrogen. 
The errors in «h this way derived are of the ord er of 50% be- 
cause of the influence of other absorption eff'ects (iDumm et al. I 
11999'). For the purpose of this work, we used values of hh for 
the above-mentioned symbiotic stars during their outbursts at the 
orbital phases between 0.3 and 0.6 to eliminate a larger influ- 
ence of the neutral wind from the giant (Table 1). Nevertheless, 
the measured values of ne can still contain a contribution from 
the giant's wind, which is in the neutral form at/around the 
orbital plane during active phases, because the ionizing pho- 
tons are blocked by the neutral disk-like formation suiTound- 
ing the hot star. For the mass loss rate from the giant of a few 
xlO"^ M0yr ' (ISkopalll2005l) . the column density of its wind 
in the direction to the hot com ponent at the orbital phase 0.5 and 
/ = 90° is of ^ 1 X 10^2 cm-2 jsekeras & Skopal 2012b). 



In spite of large errors in hh and R^, often derived from non- 
simultaneous observations. Fig. |4| shows that the nu values are 
consistent with those calculated within the WCZ model for pa- 
rameters of the hot star and its wind derived from observations 
during active phases. 
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Fig. 5. Comparison of the model emission measure calculated 
according to Eq. ( 1211 1 throughout the Hn region as a func- 
tion of the radial distance r (solid line) with the values derived 
from the spectra of selected symbiotic stars during active phases 
(Table 1). 9(H ii) is the opening angle of the H n zone. The same 
model as in Fig. |4] was selected. 

3.3. Emission measure of the ionized wind from the hot star 

In this section we calculate the emission measure, EM, of the 
ionized part of the hot star wind from regions above/below the 
neutral disk-like zone and compare it to that derived from obser- 
vations during active phases. 

Modelling the SED of symbiotic binaries during active 
phases showed that the emission measure of the nebular com- 
ponent of radiation i s between a f ew xlO^^ and a few xlO*" 
cm -3 (see Table 4 of lSkopall llOO?). In addition, the modelling 
revealed the presence of a low-temperature nebula (LTN, ~ 
14 000 K), which is subject to eclipses, and a high-temperature 
nebula (HTN, Tg > 30 000 K), which is seen during eclipses as 
the only component. Later, Skopal (2006) ascribed the LTN to 
the ionized wind from the hot star, which develops during active 
phases. In our model, EM is determined as 



EM ^ \ A^^(r,6»)dy, 
Jv 



(20) 



where A^h('", O) is the density of the ionized hydrogen given by 
Eq. ([TJ, calculated throughout the volume V of the Hii zone. 
In the spherical coordinates (r, 6, 0), the volume element dV - 
sin6drd0d(p, and in the WCZ model, the density distribution 
is azimuthally symmetric (i.e. along the coordinate 0). This al- 
lows us to express Eq. ( l20l l as 



EM = 471 



/- oo /- e(H II) 
Jr. Jo 



N^(r,e)r^ sin 6 drde. 



(21) 



where 0(li ii) is the 6 coordinate of the H i/H ii boundary at the 
given r. Figure |5] shows the calculated EM of the ionized part 
of the wind as a function of r for the same model that we used 
to calculate rtu in Sect. 3.2. It shows that the observed values 
of the emission mea sure, as derive d by modelling the SED (i.e. 
£Mltn in Table 4 of ISkopal II200I Table 1 here), are in a good 
agreement with the model prediction. We excluded the case of 
BE Cyg, because it was not possible to extract the LTN compo- 
nent from the SED. A problem in comparing the observed and 
modelled emission measure is connected with a shielding effect 



Fig. 6. Model of the ionization structure of the hot components 
in symbiotic binaries during active phases (Sect. 2, Eq. ( fT3] )) 
as seen on a cut perpendicular to the orbital ix,y) plane con- 
taining the WD, which spins along the z-axis. The neutral zone 
has the form of a flared disk. It is optically thick to the distance 
Rii, i.e. the radius of the warm pseudophotosphere. The ionized 
zone has the opening angle ^(Hii). The model corresponds to 
Rf ~ II Rq (see Eq. The inplot figures the central WD 

as a black circle with the radius /?wd, which is the source of 
the optically thick wind. The wind becomes optically thin at the 
distance (the shadow circle), which simulates the WD's pseu- 
dophotosphere, the source of ionizing photons. 



for objects with a high orbital inclination. In the real case, it 
is probable that the outer rim of the flared disk will shield the 
inner parts of the ionized wind. As a result, the calculated emis- 
sion measure EM^ai > /sMltn- However, we did not include the 
shielding effect in the calculation, because the orbital inclination 
and the opening angle of the neutral zone are not well-known 
parameters. 

According to Fig.|5] we can conclude that the emission mea- 
sure of the ionized wind from the hot star, calculated within the 
WCZ model (Eq. (ITTTi). is consistent with that derived from ob- 
servations by modelling the SED of symbiotic binaries during 
active phases. 

4. Conclusion 

We investigated the ionization structure around the hot stars in 
symbiotic binaries during active phases. Applying the wind com- 
pression model (Sect. 2), we found that a neutral disk-like zone 
can be formed as a result of the enhanced wind from the ro- 
tating WD at its equatorial plane. Around the pole directions, 
above/below the neutral disk, the wind is diluted by a fast rota- 
tion, which makes it easily ionized by the central hot star (i.e. 
the WD's pseudophosphere in the model). Examples of the cor- 
responding Hi/Hii boundaries are shown in Fig.[T] A represen- 
tative model of the ionization structure of the hot component in 
a symbiotic binary during active phase is shown in Fig.|6] 

Applying the wind compression model to physical param- 
eters of the hot stars in symbiotic binaries, we found that the 
neutral disk-like zone can be created during active phases, when 
the mass loss rate enhances to ~ 2 x lO^^Mgyr"'. However, 
during quiescent phases, the formation of the neutral region is 
unlikely because of insufficient mass loss rate only of a few 
xlO-** Moyr-i (Sect. 3.1., Fig.|2]i. 
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We verified the model by comparing the hydrogen col- 
umn density, calculated in the radial direction from the hot star 
throughout the Hi zone and the emission measure of the Hn 
zone, with quantities derived from observations (Sects. 3.2. and 
3.3.). Figs.|4]and|5]show that both parameters are consistent with 
those calculated according to the wind compression model for 
M = 2 X lO^^Moyr"' and = 2 000km s ', which represent 
mean values derived from observations during active phases. 

In view of these results, we propose that the flared optically 
thick disks, indicated in the spectrum of active symbiotic bina- 
ries with a high orbital inclination, can form from the enhanced 
wind, which is compressed towards the equatorial plane due to a 
fast rotation of its source, the central WD. In this way we justi- 
fied the ionization structure of the hot components in symbiotic 
binaries during active phases, as we had suggested in pre vious 
papers (e.g. ^kopal 2005, 2006; Skopal et al. 2006. 1201 ih . 
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